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مقایسه‌ی روش هم‌هندسه با اجزای محدود در تحلیل دینامیکی مسائل کشان دوبعدی * 


بهروز یدزی 7 احسان ژیانی عید گاهی ٩(‏ علیرضا حسن‌زاده طاهری(۳ ناصر ظریف مقدم ا؟ 


۲ 
چکیده مرضوع این مقاله تحلیل دینامیکی مسائل دوبعدی کسان به روش هم هندسه و مقایسه یآن با روش اجزای محدود است. برای 
این منظور رابطه‌سازی معادله‌ی دیفرانسیل حاکم با استفاده از توابع پایه اسپلاین‌ها به‌دس ت آمده است. با طرح چند مثال که د رآن از 
ماتریس‌های جرم سا زگار و متم رکز استفاده شده است, به مقایسه‌ی جواب‌های روش‌های اجزای محلود و هم‌هناسه پردانته شلده است. 
همچنین تأثیر افزایش تعداد نقاط وارسی و افزایش درجه منحنی اسپلاین‌ها بر دفت جواب مورد بررسی قرا رگرفته است. به علاوه. جون در 
تحلیل دینامیک ی گذرا از روش نیومارک استفاده شدهء به تأثیر عامل‌های 6 و بر حل مساله پرداخته شده است. نتایج حاصل, دقت بالای 
روش هم‌هندسه را در مقایسه با روش اجزای محدود. علی‌رغم بهکارگیری تعداد بسیا رکم‌تری از درجات آزادی و در نتیجه کاهش ابعاد 

دستگاه معادلاات و زمان حل مسأله» تال می‌دهنك. 


واژه‌های کلیدی روش هم‌هندسه تحلیل دینامیکی گذرا. اسپلاین‌ها؛ ماتریس جرم سازگار و متمرکز. 
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* تاریخ تصویب مقاله ۸٩/۰۵/۱۷‏ و تاریخ دریافت نسخه نهایی اصلاح شده ٩۱/1/۱۳‏ 
(۱) نویسنده مسئول: دانشیان گروه مکانیک دانشگاه فردوسی مشهد. 

(۲) دانشجوی کارشناسی ارشد سازه. دانشگاه صنعتی شاهرود 

(۳) دانشجوی کارشناسی ارشد. مکانیک دانشگاه فردوسی مشهد 

(4) دکترای سازه دانشکده فنی شهید منتظری مشهد. 


مقد مه 

امروزه تحلیل مسائلی که فرا روی مهندسان و محفقان 
قرار می‌گیرد. آنان را ملزم به حل معادلاتی می‌کند که 
اغلب به شکل دیفرانسیلی هستند. در بیشتر مواقع حل 
این معادلات به روش‌های تحلیلی مقدور نیست و 
به‌ناچار باید از روش‌های عددی از قبیل روش 
تفاضل‌های محدود (۱۷600 1۱111676006 ۳10166) 
روش اجزای محدود (۱۷۲6/۵۵۵ )حعصعاظ فاص 
روش‌های بدون جزء (16/000 عوعلهع]۷)؛ مانند 
روش نقاط محدود (۱6۵0 تصلوظ عالطا و 
روش‌های جدیدتر؛ مانند روش تحلیل هم‌هندسه 
(۱۷۲۵00 ماو رامش عنتاهط60ع190) استفاده کرد. این 
روش‌ها هر یک در پی دیگری آمده و دو هدف اصلی 
را دنبال کرده‌اند. این اهداف عبارتند از اولاً ایجاد 
کارایی» دقت و سرعت بالاتر و انیا ایجاد امکاناتی 
جدیدتر برای حل مسائل و رفع مشکلات روش‌های 

استفاده از روش اجزای محدود در تحلیل 
سازه‌های پیچیده و مکانیک محاسباتی به امری معمول 
تبدیل شده است. شاید بتوان گفت که توسعه‌ی این 
روش در سال‌های گذشته و وجود نرم‌افزارهای متعدد 
تجاری سبب کاربرد وسیع این روش در علوم مختلف 
هقی سازه هن مکانک مایا شاه اسر اس 
باید دانست مزایایی مهم مانند نظام‌دار بودن روش 
اجزای محدود در یافتن توابع تقریب‌زننده» مدل‌سازی 
شکل‌های پیچیده. قابلیت تبدیل به الگوریتم‌ها و 
برنامه‌های رایانه‌ای و پیشرفت سریع رایانه‌ها؛ از علل 
فراگیر شدن استفاده از این روش می‌باشد. اما روش 
اجزای محدود دارای معایبی نیز هست که از آن جمله 
می‌توان به عدم تطابق کامل مدل تحلیل( یعنی شبکه 
اجزای محدود) با هندسه‌ی موردنظر و خطای ناشی از 
آن اشاره کرد. همچنین صرف هزینه و وقت قابل 
توجه برای تولید شبکه‌بندی مناسب از جمله 


مشکلات روش اجزای محدود محسوبت می‌شود. 
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مقایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی... 


گفته می‌شود که به طور میانگین برای حل یک مسأله 
اجزای محدود حدود هشتاد درصد از زمان حل ی 

روش عددی حل مسائل با استفاده از اسپلاین‌ها 
روشی جدید است که مفاهیم آن در سال ۲۰۰۵ برای 
اولین بار توسط هیوز(65ع۳10 1.1.5) و همکارانش 
با نام روش هم‌هندسه معرفی گردید [1]. این روش 
واه ایس یگدز که هلسسین او هلال تز لیس 
جواب‌های آن با استفاده از اسپلاین‌ها و نوع توانمندتر 
آن یعنی تربز (61106-ظ لقصمتای؟ حصصمئتصنا-صمل) 
می‌باشد. ایده‌ی اصلی این روش بر استفاده از همان 
توابع پایه نربز مورد استفاده برای تولید هندسه وه 
منظور تقریب پاسخ معادله‌ی دیفرانسیل استوار است. 
به اين ترتیب می‌توان گفت که با ابداع این روش 
فرایندهای طراحی و تحلیل در یکدیگر ادغام شدند 
که این امر می‌تواند پیشبرد عظیمی در زمینه‌ی 
روش‌های عددی به شمار آید. تا کنون این روش در 
حوزه‌های مختلف مکانیک محاسباتی؛ از جمله تحلیل 
سازه‌ای» اندرکنش سبال-سازه. انتشار حرارت. 
ارتعاشات سازه‌ها و غیره مورد استفاده قرار گرفته 
است [درم2]. در ستال ۳۰۰۹ این روش برای تحلیل 
ارتعانانت. ازاد‌سانه‌های ماو ماه شمان ها 
تیرهاء صفحه‌ها و پوسته‌ها مورد استفاده قرار گرفت 
[4]. در مقاله‌ی نام‌پرده با مقایسه‌ی نتایج حاصل از 
روش ایزوژئومتریک و اجزای محدود. با استفاده از 
تعداد یکسانی از درجات آزادی. دقت و توانمندی 
ایین روش در مقایسه با روش اجزای محدود در 
تحلیل ارتعاشات آزاد سازه‌های مختلف» نشان داده 
شده اسست. همچنین در سال ۲۱ روش 
ایزوژئومتریک توسط ابوالبشری و همکاران [1] برای 
تحلیل ارتعاشات آزاد مسائل تنش صفحه‌ای به کار 
گرفته شد. در این پزوهن برای, نجستین بار دفست ق 
کارآمدی روش هم‌هندسه در مسائل وابسته به زمان 
مورد بررسی قرار گرفته است. 

با دیدگاهی قدری متفاوت. و در واقع کلی‌تره 
نگارندگان این مقاله روشی را دنبال کرده‌اند که مستقل 


سال بیست و چهارم» شماره یک» ۱۳۹1 


بهروز حسنی- احسان ژیانی- علیرضا حسن زاده- ناصر ظریف 


از هیوز و همکارانش بوده است و علاوه بر برخی 
مسائل معمول. در حل معادلات دیفرانسیلی که 
ضرایب دیفرانسیل‌ها خود توابعی هستند. با موفقیت 
به کار گرفته شده است [۸.و7]. در این دیدگاه به 
جای اندیشه‌ی جایگزین کردن جزء‌های اجزای 
محدود با جزهء‌های گرهی(61676005 16006 حل 
مسأله برای هریک از مولفه‌های میدان به صورت یک 
سطح (يا شبه سطح در ابعاد بالاتر از دو) تصور شده 
است و مختصات نقاط وارسی آن, به عنوان مجهولات 
مسأله در نظر گرفته شده‌اند. در اين مقاله به‌کارگیری 
این روش و کیفیت حل به دست آمده در مسائل 
دینامیکی مورد بررسی قرار گرفته است. 

برای انجام تحلیل گذرای دینامیکی؛ در این روش 
نیز مشابه روش اجزای محدود. ماتریس‌های جرم. 
سختی و میرایی (در صورت وجود) ساخته شده 
است؛ هر چند که ماتریس‌های فوق مربوط به نقاط 
وارسی می‌باشند و فاقد مفهوم فیزیکی هستند. برای 
تابع اولیه‌گیری» زمانی از روش نیومارک که یک روش 
تابع اولیه گیری مستقیم است استفاده می شده است. 

در بخش دوم این مقاله به معرفی تحلیل 
دینامیکی و در بخش سوم به معرفی اسپلاین‌ها و فن 
نربز و شرح مختصر عامل‌های آنها پرداخته شده 
است. بخش چهارم روش هم‌هندسه را در مقایسه با 
روش اجزای محدود معرفی کرده و در بخش پنجم 
نحوه‌ی به دست آوردن رابطه‌سازی معادلات حاکم. 
به روش اسپلاین‌ها ارائه شده است. در بحش ششم با 
ال اتکی )تانب و تا وم فان وی تون 
در بخش هفتم نتایج و پیشنهاداتی در این خصوضن 
بیان شده است. 

تحلیل دینامیکی 

معادله‌ی حرکت سازه در نوسان اجباری در اثر اعمال 
نیرو پس از انجام گسسته‌سازی دامنه مسأله به صورت 
زیر است که در آن 16 ماتریس جرم. لا بردار شتاب؛ ) 
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لآیردار جابه‌جا باشد [9]. 
ی وا بو ی 


1 + 00 + 1 < ۴ 0) 


دستگاه معادلات فوق از لحاظ ریاضی بیانگر یک 
سیستم معادلات دیفرانسیل خطی مرتبه دوم است که 
برای حل آن می‌توان از روش‌های جمع آثار مدها و 
تابع اولیه‌گیری مستقيم استفاده کرد. در اینجا روش 
دوم مورد نظر بوده است. برای حل این معادلات 
می‌توان از روش‌های مختلف عددی؛ نظیر روش‌های 
اجزای محدود. تفاضل‌های محدود و یا روش‌های 
بدون شبکه به منظور گسسته‌سازی دامنه و نیز 
تفریب‌سازی و درولن‌یابی بهره جست. در همه‌ی این 
روش‌ها برای داشتن دقت کافی به کارگیری نقاط 
گسسته‌سازی نسبتاً زیادی مورد نیاز بوده و در نتیجه 
ماتریس‌های ضرایب به دست آمده بزرگ و پر هزینه 
می‌باشند. نظر به این ویژگی روش هم‌هندسه که در 
آن امکان تعریف سطوح با نقاط وارسی محدود میسر 
می‌باشد در اینجا این روش مورد استفاده قرار گرفته 
است. 

برای تابع اولیه‌گیری زمانی از یک روش عددی 
گام به گام با استفاده از خانواده‌ی روش‌های نیومارک 
بهره گرفته‌ايم که خلاصه‌ی آن در جدول (۱) قاببل 
مشاهده است [9]. 

در این جدول. » و ۷ عامل‌هایی هستند که می 
توان آنها را برای به دست آوردن دقت تابع اولیه‌گیری 
و پایداری روش تعیین کرد. تاکنون روش‌های مختلفی 
قفاب اتقات ای گوتا کرت یا رارسا 
پيشنهاد شده‌اند و هر یک از آنها مزایا و معایب 
متقافی داش سوارو اراشه هه دی ورن و 
پرکاربردترین روش‌های پیشنهادی مورد استفاده در 
تحلیل‌های گذرای دینامیکی به شمار می روند که از 
پایداری و دقت مناسبی برخوردار هستند. از این‌رو در 
هه فا ند از هی تسا افاده نسته این 
جزئیات بیشتر در مورد پایداری و دقت هر یک از این 


روش‌ها در مرجع [9 آمده است. 
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جدول ۱ برحی از روش‌های خانواده‌ی نیومارک 


روش شتاب میانگین 
4 2 2 26 < ۲۷ 2 < 0 
روش تفاضل مرکزی | 8/5 2 26 < ۷ | 3/2 <0 
روش تفاضل پسرو 26-2 < ۷ 2 - 0 


معرفی اسپلاین‌ها و فن نربز 
اسپلاین‌ها در واقع ابزاری ریاضی هستند که می‌توان 
با استفاده از اطلاعاتی اندک. منحنی‌ها و سطوح 
پیچیده را به صورت عددی مدل کرد [10,11]. برای 
تولیل یک متس مایت حوسطایت ای بات با 
برخی از تعاریف و مفاهیم آن آشنا بود که به احتصار 


در ادامه آمده است. 


بدا گرد ردان اساوی ات تا لا دز نظ 
گرفته می‌شود. این بردار شامل مجموعه‌ای از اعداد 
حقیقی است؛ به طوری که رابطه‌ی رررلاگ رل و 
(2,...,90-1یا,10) در آن برقرار است. بردار 7 
پردار گره (۷۰6۵۵۲ 16006) و زلاها مقادیر گرهی (000) 
نامیده می‌شوند. بازه (ربزلا مزا | را یک دهانه 
گرهی ٩020(‏ 01 می‌نامیم که می‌تواند طول آن 
صفر باشد؛ زیرا گره‌ها لزوماً دارای مقادیر متمایز 
نمی باشند. اگر فاصله‌ی بین گره‌ها مساوی باشد. در 
ایسن حالت بردار گره را یکنواخت (طه‌گنوتا) 
می‌نامند. در غیر این صورت این بردان غیر یکنواخت 
(صمئتصلا-ظ0) خواهد بود. 


توابع پایه اسپلاین. آامین تابع پایه اسپلاین با 
درجه‌ی ‏ (یا مرتبه 1+) رابا () ,رل نشان 


> > را ا 1 


)۲( نت - (۲) و 


اوت ۹ 0 
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توجه شود که برخی از خواص عامل‌های موجود در 
روابط فوق عبارتند از: 

ِ- )۳ ورلا یک تابع پله‌ای (۳۱01100 8000) بوده و 
به جز در بازه‌ی (,,بلا,]> دا مقدار آن در سایر نقاط 
صفر است. 

۲-برای 0 < و آن‌گاه () ما ترکیب خطی از دو 
تابع پایه با درجه‌ی 0-1 می‌باشد. 

۳- برای محاسبه‌ی کلیه‌ی توابع پایه به بردار گره تا 
و درجه 0 نیاز می‌باشد. 

6- در صورتی که رابطه‌ی(۱-۲) به تقسیم م9 برسده 
آن را برابر با صفر لحاظ می‌کنند. 

۵- )۲ ۳۳ ها جند جملهای‌های تکه‌ای 
(ادتصجمطوا۳0 ۳1۵۵۵۷۷156) هستند که بر روی یک 
حط حقیقی تعریف شده و فقط بازه لاملا ] را 


تحت تأثیر قرار می‌دهند. 


تولید سطوح اسپلاین و نریز. می‌توان یک سطح 
نربز را که دارای درجه‌ی 0 در جهت ۷ و درجه‌ی 
و در جهت ۷ می‌باشد را با رابطه‌ی(۳) به شکل زیر 
محاسبه کرد [11]. 


۱ 


رده زار () و (ن) «ر۵ 2 


٩),۷(- ِِ 


0۷1 ت 
ِ ۷ 1 ( ۳ "۳ م2۵ 


0ز 0-ز 


)۳( 


که (,,۳) شبکه نقاط وارسی و (ر,۳) وزن مربوط 
به آنها می‌باشند. توجه شود که اگر وزن کلیه نقاط 
مساوی یک لحاظ شود آن گاه منحنی‌ها و سطوح 
تولیدی را با عنوان اسپلاین می‌شناسند و در غیر ایین 
همچنین ((),,() و ((۷),,ل() توابع پایه غیرگویای 
اسپلاین (100اعص وتععط عصتامو-ظ اقعمتاه: -عمل) 
می‌باشند که بر روی بردارهای (۵و4) تعریف 


می‌شوند. 
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بهروز حسنی- احسان ژیانی- علیرضا حسن زاده- ناصر ظریف 


3 (اوسواو یم لا وس و ولا ,10:۰ < ل] 
1+ + 

(6) [لوسلویو۷ وی ۷و0و0 ۷ 
اجه 11 


9 ۲ د که روابط 1++(<۲ و 
1+و+و عو نیز برقرار می‌باشند و اگر (۷,) 13۳7 
به شکل رابطه‌ی(۵) در زیر تعریف شود: 

622 (0) از 
0 ۱ کب (۷.) 18 


۲۳۳ م, ‏ () معا 2 2 
اکنون می‌توان رابطه‌ی(1) را به شکل (۷ باز نویسی 
کرد. 


۹2 


با استفاده از رابطه‌ی(۷) می‌توان هر گونه سطح و 
رویه‌ی پیچیده‌ای را ایجاد کرد که در این‌جا در روش 
هم‌هندسه برای حل مسائل از این روابط برای 
الگوسازی جواب‌ها استفاده می‌شود. قابل ذکر است 
که در مدل‌سازی و تحلیل هندسه‌های پیچیده‌تر؛ از 
جمله مقاطع مخروطی, به منظور مدل‌سازی دقیسق 
هندسه بایستی از منحنی‌های نربز به جای بی 
اسپلاین‌ها استفاده کرد. در این مقاله رابطه‌سازی مسأله 
در حالت کلی و بر پایه‌ی منحنی‌های نربز استخراج 
شده است؛ ولی از آن‌جایی که مثال‌های حل‌شده 
هندسه‌ی ساده‌ای دارند. برای تشادن مسا اما از 
بی‌اسپلاین‌ها استفاده شده است. اطلاعات تکمیلی در 
این خصوص از مراجع [11,12] قابل دست‌یابی 
می‌باشد. 


روش هم منت در مفایسه با روش احزای 
محدود 
روش هم‌هندسه. در اصل از ترکیب فن‌های طراحی 
توسط رایانه (1(69120 ۸060 ععاامه0)) و روش 


اجزای محدود حاصل شده است [1,11]. چهارچوب 


سال بیست و چهارم؛ شماره یک ۱۳۹۱ 


کلی تحلیل هم‌هندسه با استفاده از فن نربز در مقایسه 
با روش اجزای محدود. مبتنی بر اصولی می‌باشد [1] 
که عبارتند از: 

۱- در این روش به جای جزءها. شبه جزء‌ها وجود 
دارد که با ضرب بردارهای گره حاصل می‌شوند و 
باید توجه کرد که به هیچ عنوان مفهوم آن مانند 
جزءها در روش اجزای محدود نمی‌باشد. به عنوان 
مثال‌در مسائل دوبعدی مانند مسائل تنش /کرنش 
مسطح این شبکه جزء به صورت 1*۷ خواهد بود 
که پیش‌تر به آن اشاره شده است. 

۲- جزء‌ه | از تقسیم دامنه توسط دهانه‌های 
گرهی(0005٩‏ 101) به دست می‌آیند. 

۳- هندسه‌ی مسأله با مشارکت شبکه نقاط وارسی و 

توابع پایه تولید می‌شود. 

۶- با استفاده از مفهوم هم‌عاملی (ماتاعصطه‌تهووو]) 
مجهولات را مانند تغییرمکان‌ها» تنش‌هاء سرعت. 
حرارت و غیره با همان توابع پایه‌ای» که برای تعریف 
هندسه استفاده شده. بیان کرده استت: 

۵- درجات آزادی عبارتند از متغیرهای وارسی 
(معاطمتت۷2 0081۲01)) که همان ضرایب توابع پایه 
قوسا خل 

در جدول (۲) می‌توان مقایسه‌ای را بین روش 
اجزای محدود و روش هم‌هندسه از دیدگاه مشابهت 
برخی مشخصه‌های کلیدی در هریک از دو روش 

ملاحظه کرد. 


به دست آوردن رابطه‌سازی روش هچ تاد سید 

در ایین بخش پس از دست‌یابی به تجربیاتی در 
چم وص روش هم‌هندسه در کارهای پژوهشی 
گذشته [۸ و7]. ابتدا معادله‌ی دیفرانسیل حاکم نو 
مساله‌های تنش/کرنش مسطح معرفی شله و سپس 
مراحل دست یابی به فرمول‌بندی روش هم‌هندسه 
برای حل این مساثئل بیان شده است. 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


جدول ۲ مقایسه‌ی روش اجزای محدود و روش هم‌هندسه 


روش اجزای محدود روش هم‌هندسه 

متغیرهای گرهی متغیرهای وارسی 

مقادیر گره‌ها در بردار گره 
یر ر بردار گر 


عدم انجام درون‌یابی نقاط و 


انجام درون‌یابی نقاط و 


متغیرهای گرهی با توابع پایه | متغیرهای وارسی با توابع پایه 


توابع پایه از نوع چندجمله‌ای توابع پایه نربز 
زیردامنه‌ها وصله‌ها 


همان‌گونه که از نگره‌ی کشسانی برمی‌آید. 

معادلات دیفرانسیل (۸) به طور هم‌زمان بر مسائل 
0 0 0۷ 0 0 

۳ ۳ ۳ 017 5 س 9 05 ده) 0 


[5-۵ 


0( 
00 0 07 00 0 
+ )2 - |( + ) مه - 
02 1 07۱ 
2 ۲۷ 2 


روابط )٩(‏ بیان می‌شوند. 


۷ 00 ۱ ۳۲ 00 ۳ 
بّ 2 2 ۱۳ 2 ۳ 2 6 


00 1 00 
10 وچ اج ۵ 4 8 22 2۵ را 


که در این رابطه 2 و 7 بردارهای یکه در جهعت 


8) 


محورهای »و ۰۷و با و با مژلفهه ای افقی و 
عمودی نیروهای گسترده سطحی وارد بر جسم 
هستند. در این روابط. نها ضرایب کشسانی 
می‌باشند که با فرض ثابت بودن ضریب پواسون و 
ضریب کثسانی در سراسر دامنه‌ی حل مسأله, در 
حالت کلی برای مسأله‌های تنش مسطح با مصالح 
همگن از روابط (۱۰) به شکل زیر 


تشریه مهندسی عمران دانشگاه فردوسی مشهد 


مفایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی... 


9 ۱۰12 ۳ < وو) < 011 


مد < 666 , 6 


و در مسأله‌های کرنش مسطح با مصالح همسان گرد 
از روابط (۱۱) به دست می‌آیند [13]. 


شم لاش - 
22 ۳ 1-۱-22 

1 )1-۶ 

۳ ۱1) 

(1-۳-22) رن 17) 


طبر 
نا ع< وم رسسسل. 
12 رورت 


11 


رادشه و اف مب له ببه:؛ هجوت ووانط (۲۲) 
برای تغییر مکان (۱۸6,5 و (۱۳) برای تغییر مکان 
(۷)۲,5 در دامنه‌ی مسأله در نظر گرفته شده است که 
در واقع از تعریف تولید رویه‌ی اسپلاین حاصل شده 
است. همچنین تغیبر مکان در (2,9) متناظر با هر نقطه 
(10,) در دامنه‌ی مسأله از روی ارتفاع رویه حل 
استخراج می‌شود که ۲ و 5 همان عامل‌هایی هستند 


که مقادیرشان با توجه به بردار گره حاصل می‌شود. 


) ۲ 0 ۷ (ه) ((2) ,۷( 2۵ 


۵ (۷, )( 9 ۳ 
2۸۵ ):9(۷, 


مز 0-ز 


همچنین در کلیه روابط بایستی به جای ۶ و ۷ 
نیز از روابط عاملی مشابه استفاده کرد و آنهارا با 
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عامل‌های ۲و جایگزین نمود. در این رابطه" و 5 
مختصات مربوط به دستگاه مختصات پارامتری 
هستند. که در روابط (۲-۷) با لاو 7 نشان داده 
شده‌اند. در این جا به دلیل نمایش مولفه‌های متغیر 
میدانی با لا و ۷ از این عامل‌ها استفاده شده است. 
اکنون می‌توان تابع‌نمای (لمدمنا۳۱0) مربوطه را به 
شکل رابطه (۱۶4) نوشت. 


9 ( 1 (ندره) 29 - ۲1 


که در اين رابطه لا بردار تغییر مکان و (3)00,۷ و 
(1)۷ به ترتیب بخش‌های دوخطی (7هعطنا3) و 
خحطی (110687) تابع‌نسا طبق روابط (۱۱و۱۵) 


2 2 
(2) ۳ یط عمط > 10 ,)۳ 


ِ نع 1-2 ()) 66 + 62 ووع + ِ م26 


(۱۵) 0106[ ا( 62۵ 


+ نام - ,۷6 + ۳3 " - ()1 
( ۷,(0۲ + ی طمعحم - 060 [(۷ 


زیرنویس 08100 در این‌جا مشابه روش اجزای 
محدود که برای المان‌ها استفاده می‌شود. برای 
وصله‌های تشکیل‌دهنده دنه مساله اسفاده»شسته 
است. همچنین از پارامتر 7 برای انتگرال‌گیری روی 
مرزها و از پارامتر # برای انتگرال گیری روی سطح 
استفاده شده است. ما و با نیروهای خارجی 
(107665 ۲201100) و پارامتر 8 نیز ضخامت هر 
وصله می‌باشد. 

پس از نوشتن شکل ضعیف معادلات و محاسبه 
ژاکوبین‌های مورد نیاز (۲) دستگاه معادلاتی به شکل 
رابطه (۱۷) به دست می‌آید که با حل آن می‌توان 


مختصات ارتفاع نقاط وارسی را برای تولید رویه‌ی 


سال بیست و چهارم؛ شماره یک ۱۳۹۱ 


جواب محاسبه کرد. سپس با استفاده از رابطه تولید 
برای هر نقطه دلخواه از مسأله و یا ترازهای تنش و 


تشی کان زا فراع شا له نه مت ارو 


][)1۱ + ][)0( < )۳( ۱۷ 


که در این رابطه[10] ماتریس جرم (ل) بردار 
قانت هام هش ی اف ۴ ردان مها و 
فاشتکا نها امه قایل خهتر سیخ کته 
مقادیر نیروهای معادل بر روی نقاط وارسی ناشی از 
اعمال بارهای گسترده‌ی سطحی و حجمی به جسم 
با استفاده از جمله‌های رابطه‌ی تابع اولیه‌گیری (۱) 
به‌سادگی قابل محاسبه است. بردار (۳) در رابطه‌ی 
(۱۷) بیانگر مجموع این نیروها است. قابل ذکر است 
که هر چند در حالت کلی لزوماً نقاط وارسی بر روی 
هندسه‌ی جسم قرار ندارند. ولی این امر به معنای 
اعمال بارگذاری و پا شرایط مرزی تکیه‌گاهی به 
خارج از فضای هندسی نیست؛ و همان‌طور که اشاره 
شد. برای اعمال شرایط مرزی موجود تنها کافی است 
سهم هر یک از نقاط وارسی با استفاده از رابطه‌ی 
(۱1) محاسبه شود. همچنین سهم هر یک از نقاط 
وارسی در نتیجه‌ی اعمال بارهای متمرکز نیز با استفاده 
از همین رابطه و به کمک مفهوم تابع دلتای دیراک 
(ممناعصیط ماه1۲ 20تذ) به‌سادگی قابل محاسبه 
است. اعمال شرایط مرزی ضروری همگن نیز مشابه 
روش اجزای محدود و بر روی نقاط وارسی مرزی 
صورت می‌گیرد. اضافه می‌ گردد که تنها اعمال شرایط 
مرزی ضروری غیرهمگن به دلیل عدم ارضای 
خاصیت تابع دلتای کرونیکر (1۲(6۱62 عمعامع08ک1) 
توسط توابع پایه نربز. در این روش برخلاف روش 
اجزای محدود به‌صورت دقیق امکان پذیر نیست. 
ماتریس جرمی به دست آمده از رابطه‌ی (۱۷) 
ماتریس جرمی سازگار نامیده می‌شود که در حالت 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


کی شک دزی کر انست: دز میبایل زا تاه 
درجات آزادی زیاد. به منظور افزایش سرعت 
محاسبات. ماتریس فوق را با استفاده از روش‌های 
ملق قظرخ مس کفی. تکن از انم ووتن‌ها که در این 
مقاله نیز مورد استفاده قرار گرفته است. قرار دادن 
مجموع عناصر روی هر سطر به جای عنصر نظیر 
روی قطر و صفرکردن ساير عناصر روی سطر است 
اتف داش کی و توافمته و تفت باتش 
در روش مذکور با حل نمونه‌هایی نشان داده شده 


تست 


حل جند مثال 
به واسطه‌ی عدم وجود مثال با حل دقیق در کتب 
علمی در این‌جا مثال‌هایی ارائه شده که در آن به 
مقایسه‌ی حل به دو روش اجزای محدود و هم‌هندسه 


پرداخته شده است. 


مثال۱. صفحه طره‌ای شکل زیر. که تحت اثر نیروی 
هارمونیک برشی 109) 200605 < (۳)0 در انتها آزاد 
می‌باشد در نظر گرفته شده است. طول صفحه ۱۰ 
سانتیمتر و عرض و ضخامت آن به ترتیب ۲و ۱ 
باتتیجفر غی پاش تخیر مکتان وانسکه به زشان: که 


۶ 


۱( 
بر 


۹۹۹۹۹۹ 
۱ ۵ 
4 0 ) 
۳ 
لح ۲-20 


۷ / 


1000۳ 


شکل ۱ شبکه‌بندی اجزای محدود مثال حل‌شده۱ 


تحت 


خواص صفحه‌ی فوق به شرح زیر است: 


ِ ۳ 18 ِ 
< م , 0.15 < ۷ 9 0 << ۲ 


5 -< ۸ , خت 1 
ار 


3 
مسأله‌ی فوق با روش هم‌هندسه و با تعداد ۵۵ 
نقطه‌ی وارسی مطابق شکل (۲) می‌باشد. البته دقت 


تشریه مهندسی عمرال دانشگاه فردوسی مشهد 


مفایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی.. 


شود که در شکل (۲) هر مستطیل یک جزء نیست و 
اصلاً در این خصوص شباهتی بین اجزای محدود و 
روش هم‌هندسه وجود ندارد. در واقع جزهء‌های روش 
هم‌هندسه از حاصل‌ضرب بردارهای گرهی به دست 
می‌آید. در حالت کلی گره‌ها در اجزای محدود بر 
روی دامنه قرار می‌گیرند در حالی که در روش 
هم‌هندسه تقاط وارسی می‌توانند بر روی دامنه نباشند. 

بردار گکره برای جهت طولی برابر با 
[او او او7۵ »و۵۰ .و۲۵ »و »و عو»] و برای جهت 
عرضی برابر با (اواواوات و۰۳۳ وتو موه) در 
نظر گرفته شده است. همچنین درجه‌ی توابع پایه 
اسپلاین در این مسأله مساوی دو لحاظ گردیده است. 


نتیجه تحلیل هم‌هندسه و اجزای محدود با 
شرایط ذکرشده و با ماتریس جرم سازگار و متمرکز 
در شکل‌های ( و۳) نشان داده شده است. 


(۵2۵ ۱ #مانعمه6) 5 ۵9۳0و 5 ۵۱69۳05 ۷۵ 6 ۵۲۱۸۸ 20۵۳۵ سس 2 


۷۷۵,۵0, 


شکل ۳ نمودار تغییرمکان- زمان مثال ۱- جرم سازگار 
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۷۵ و 5 هدن 5 5 واه جات ۴6۸ انشا 60۳۵00 2 


0۷ )۹۵,۵,۵ 


شکل ۶ نمودار تغییرمکان- زمان مثال ۱- جرم متمرکز 


ملاحظه می‌شود که در هر دو حالت ماتریس 
جرم جواب‌های مناسبی به دست می‌آید و عدم انطباق 
کامل نمودارها قابل قبول است زیرا با توجه به این که 
در روند حل دینامیکی از تغییر مکان. سرعت و شتاب 
در گام زمانی قبل استفاده می‌گردد. این عدم انطباق 
قابل پیش‌بینی است: از طرفی این عدم انطباق در 
بررسی این صفحه با شبکه‌بندی‌های مختلف به کمک 
روش اجزای محدود نیز ملاحظه می‌شود که در شکل 


زیر نشان داده شده است. 
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شکل ۵ مقایسه‌ی نمودارهای تغییرمکان- زمان با شبکه‌بندی‌های 


مختلف به روش اجزای محدود 


سال بیست و چهارم. شماره یکه ۱۳۹۱ 


قرارگرفته و با روش اجزای محدود مقایسه شده 


۷۸ ۱۰ 4اشا 
اع۴ ۰ 


۱ 
وا 12۸/82۰ 
۱۹ (۱۸۸)۴<2 ۳ 3 


(هعوونامه۱۱, ۱۸/۷/۴<۵ 
(55-9لا۱۱6۸, ۱۸۸۸/۴<3 


0۷)40۵,0,۵ 


شکل 1 نمودار تغییر مکان- زمان منال ۲ 


نسبت به 1 < ] جواب‌ها انطباق بهتری با روش 
اجزای محدود دارند. همچنین افزایش درجه‌ی منحنی 
به مقدار ۳ تأثیر چندانی بر بهبود جواب نخواهد 
داشت. همچنین بررسی انجام‌شده نشان داده که 
به هنگام استفاده از منحنی اسپلاین درجه ۳ تنها باعث 
بهتر شدن کیفیت نمودار تغییرمکان - زمان در نقاط 


اوج نمودار می‌شود. 


مثال۳. اثر تعداد نقاط وارسی بر دقت جواب. در این 
قسمت مثال ۱ ثابت نگه داشتن درجه‌ی منحنی 
اسپلاین و افزايش نقاط وارسی در دو حالت جرم 
متمرکز و سازگار مورد بررسی قرار گرفته که در زیر 
مشاهده می‌شود. 

همان طور که ملاحظه می‌شود افزایش نقاط 


وارسی باعث بهبود جواب خواهد شد. البته این روند 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


۱۰ 


ادامه توضیح داده شده است. 


عاصاه ۵ امنوهی اه ۲ماصت ۱ گن عمعوودی 8 


۷ )1۵,۵,٩( 


شکل ۷ نمودار مثال ۳با تعداد نقاط وارسی متفاوت و با 


ماتریس جرم سازگار 


وود معمصسا) کات ۳ 20۲۷۲0 آن ۲ع ۱۷ ۵۷ 600۲۳۵۵۷۲ 10 


0۷)40,0,( 


شکل ۸ نمودار متأل۳ با تعداد تقاط وارسی متفاوت و 


با ماتریس جرم متمرکز 


به هنگام استفاده از ماتریس جرم متمرکز اگر 
نحوه‌ی افزایش نقاط به گونه‌ای باشد که نسبت 
این روند به بهبود جواب کمک می‌کند ولی در غیراین 
صورت چنین نخواهد بود که در شکل (۱۰) ملاحظه 
می‌کنید. 


نشریه مهندسی عمرال دانشگاه فردوسی مشهد 


مقایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی.. 


مشاهده می‌شود در حالتی که از ماتریس جرم 
تنهایی منجر به جواب دقیق‌تر نمی‌شود. بلکه این 
فزایش نقاط باید به نحوی باشد که نسبت فاصه‌ی 
فزایش تقاط وارسی به ۱ دیده می‌شود که جواب‌ها 


ز مقدار صحیح خود دور می‌شوند. 


‌ ۱ 


0۷), 


شکل ٩‏ نمودار مثال۳ با تعداد نقاط وارسی متفاوت و با مائریس 


جرم سازگار 


(که ۱۸۸ 4ومسا) عاماه۳ اوتاانج که ۲وحاات ۱ 0 ععص ان ۷ انشا ۱ 8 


۱۷ )0,۵,( 


شکل ۱۰ نمودار مثال۳ با تعداد نقاط وارسی عتفاوت و با 


ماتریس جرم متمرکز 


سال بیست و چهارم. شماره یک (۱۳۹ 


به روز حسنی- احسان ژیانی- علیرضا حسن زاده- ناصر ظریف 


قو الم وزنت افایش قاط وارشی سعی باه 
است که تغییرات نسبت فاصله‌های نقاط وارسی در 
دو راستای عمود بر هم نیز مورد ارزیابی قرار گیرد تا 
از افزايش بی‌مورد نقاط که تأثیر منفضی بر جواب‌ها 
می‌گذارد. جلوگیری شود. 


قرار می‌گیرد و دو روش هم‌هندسه و اجزای محدود 
در دو حالت جرم متمرکز و سازگار با هم مقایسه 


می سود. 
اوه ۱۷ 2۵۳2۳۸) ۴6۱ ۷۷۰ ۸۱ 
(6300820.5 6 ۸۵20.5) ۱۸ سس 3 
(6۵۳۵۶0.5 6 ۸۷۵2۵.5) ۴6۱۸ 
(6۵00۵21.6 6 ۸۳21.5 ۱۵ سم و2 
(6۵00821.6 ۶ ۸۷۵21,5) ۴6۱۸۷ 
(6200۵22 ۶ ۸۷۵21.5) ۱۸۱۸ ۳ 
(6۵60822 8 ۵۸۷۵21.5) ۴6۸0 سس ۳ 
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شکل 11 مقایسه‌ی دو روش هم‌هندسه و اجزای محدود با خود 
و با یکدیگر به کمک عامل‌های آلفا و گاما - جرم سازگار 


(عهه ۱۵ موی ۴2۱۸ .۷۵۴ مشش 


(6۵۳۱۵20.5 6 ۸۱/۵۸20.5) /اها ۳ 


ط 


0۷)1۵,۵,( 


شکل ۱۲ مقایسه دو روش هم‌هندسه و اجزای محدود با خود و 
با یکدیگر به کمک عامل‌های آلفا و گاما - جرم متمرکز مثال ؛ 


سال بیست و چهارم. شماره یکه ۱۳۹۱ 


۱۲۱ 


در این مثال روش‌های شتاب میانگین ثابت؛ 
تفاضل مرکزی و تفاضل پسرو که از خانواده‌ی 
نیومارک می‌باشند با یکدیگر مقایسه می‌شوند. ملاحظه 
می‌شود که جواب‌ها در اين سه روش تقریباً به یک 
جواب همگرا می‌شوند و از طرفی دو روش هم‌هندسه 


و اجزای محدود نیز همگرایی خوبی با یکدیگر دارند. 


مثال ۵. تیر یک‌سرگیردار-یک‌سرآزاد شکل زير را در 
نظر بگیرید. این تیر تحت بار گسترده خطی در 
سرتاسر طول خود می‌باشد. طول تیر ۲۰۰ میلی‌متره 
عرض آن ۲۰ میلیمتر و ضخامت آن برابر ۵ میلی‌متر 


می‌باشد. خواص مکانیکی تیر به شرح زیراند: 


تص/ع1 7806 م , 029 > ۷, 06۳۵ 210 < ]1 


هندسه و شرایط مرزی تیر در شکل (۱۳) نشان 
داده شده است. همان‌طور که در این شکل نشان داده 
شده است» یک بار گسترده‌ی یکنواخت بر روی 
صفحه‌ی بالایی تير اعمال شده است که به صورت 
یی ماو کف مکی کر این ساله با شیر 
واسعه رسای کتره شآ ابا کی موز 


بررسی قرار گرفته است. 


م0 رو ع )0 


20 


مک ۴ میس ری موق کبک گرواز 


مثال فوق با استفاده از هر دو روش اجزای 
محدود و ایزوژئومتریک مورد بررسی قرار گرفته 
است. به منظور تحلیل مسأله با استفاده از روش 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


۱۲ 


یزوژئومترک» از یک شبکه‌ی ۲۰ در ۲ از نقاط 
کنترلی که شامل ۱۲۰ درجه آزادی است. استفاده شده 
ست. همچنین توابع پایه مورد استفاده از درجه‌ی ۲ و 
بردار گرهی به صورت یکنواخت در نظر گرفته شده 
ست. در حالی که برای تحلیل با استفاده از روض 


جزای محدود از ۶۰ المان درجه‌ی ۲ (المان 


سرندیپیتی) که مجموعا شامل ۳۳۰ درجه‌ی آزادی 
می‌شود. استفاده شده است. نمودار جابه جایی گره 
میانی انتهای تير بر حسب زمان در شکل (۱۶) نشان 
داده شده است. همان طور که مشاهده می‌شود علی 
رغم آن که در روش ایزوژئومتریک از تعداد بسیار 
کمتری از درجات آزادی استفاده شده است. نتایج 
کاملاً بر یکذیگر منطبق شده‌اند. از آنجایی که حل این 
مثال با استفاده از روش ایزوژئومتریک منجر به حل 
دستگاه معادلات بسیار کوچک‌تری نسبت به روش 
اجزای محدود می‌شود. لذا قاعدتاً زمان محاسباتی 


مورد استفاده در این روش بسیار کم‌تر خواهد بود. 
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و) ۲۳6 
شکل ۱۶ نمودار تغییر مکان- زمان مثال ۵- جرم سازگار 


و ۱ 


انتها برابر ۳ سانتی‌متر می‌باشد. ضخامت صفحه نیز 


تشریه مهندسی عمران دانشگاه فردوسی مشهد 


مقایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی... 


ترایز ۱۰ سانش مر انیت یاف منکیم هه ره 


شرح زیراند: 


دمم , ۷۶03 , 706۴2- 8 
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همان طور که در شکل (۱۵) نشان داده شده 
است» این صفحه تحت اثر دو بار نقطه‌ای سینوسی در 
نقاط کوگیا ام تشر «اواد خود که به صورت 
209) طزی < 2000 < )و1 < () با زمان تغییر 
می‌کنند. می‌باشد. شبکه‌بندی اجزای محدود و نیز 
شبکه‌ی نقاط کنترلی مورد استفاده برای حل این مثال 
به ترتیب در شکل‌های (۱۷ و *۱) نشان داده شده‌اند. 
همان طور که مشاهده می‌شود. شبکه‌بندی اجزای 
محدود مورد استفاده شامل ۸۵ نقطه‌ی گرهی (معادل 
۰ درجه‌ی آزادی) و شبکه‌ی نقاط کنترلی شامل ۳۳ 
نقطه‌ی کنترلی ( 11 درجه‌ی آزادی) می‌شود. در 
این‌جا میزان جابه جایی گره میانی انتهای تير بر 
حسب زمان مورد بررسی قرار گرفته و در شکل (۱۸) 
در بازه‌ی صفر تا پنج ثانیه رسم شده است. 

قابل ذکر است که در هر دو روش از تعداد 
یکسان ؛ نقطه‌ی گوسی در هر المان یا دهانه‌ی گرهی 
برای تابع اولیه‌گیری استفاده شده است و همان طور 
که شکل (۱۸) نشان می‌دهد. نتایج انطباق نسبتاً خوبی 
به ویژه در زمان‌های ابتدایی با یکدیگر دارند. 


سال بیست و چهارم» شماره یک» ۱۳۹1 


به روز حسنی- احسان ژیانی- علیرضا حسن زاده اصر ظر یف 


سس 


شکل ۱۷ شبکه ی نقاط وارسی منال" 
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شکل ۱۸ نمودار جابه جایی گره میانی انتهای تیر 
بر حسب زمان 
این در حالی است که ابعاد دستگاه معادلات 
حاصل از روش ایزوژئومتریک تقریباً یک سوم روش 
اجزای محدود می‌باشد و قاعدتاً منجر به کاهش 


سال بیست و چهارم؛ شماره یک ۱۳۹۱ 


۳۳ 
چشم‌گیر زمان محاسبات خواهد شد. 


نتبجه گیری و پیشنهادات 

پس از بررسی‌های انجام‌شده نتایجی به شرح زیر به 
دسا آمانن 

در حل مسائل هم ماتریس جرم سازگار و هم 
ماترسیبجرم متمر کر تجواب‌های: میتی ازابه کردند: 

ریش دزی« تین استی دینش رن 
تعداد نقاط گوسی باعث بهبود جواب در تحلیل 
دینامیکی خواهد شد. 

در تحلیل دینامیکی به کمک ماتریس جرم 
تیک افتاسی تعلاد قاط واسی بان ار بط 
همراه باشد تا به جوابی مناسب‌تر نائل گردد. و شرط 
پادشده آن است که نسبت فاصله‌ی طولی و 
عرضی(۲۵0 /۸۵060) نقاط وارسی بیشتر از ۶ نباشد. 

با توجه به نمودارهای ارائه شده در مقایسه‌ی 
روش‌های شتاب میانگین ابست. تفاضل مرکزی و 
تفاضل پسرو با یکدیگر به کمک روش هم‌هندسه 
جواب‌های مناسبی ارائه. شد. علاوه بر این؛ تحلیل 
دینامیکی با این سه روش به کمک روش‌های اجزای 
محدود و هم‌هندسه با یکدیگر مقایسه شد که 
جواب‌ها تقریباً بر یکدیگر منطبق شدند. 

عدم انطباق کامل جواب‌ها در تحلیل دینامیکی به 
کمک دو روش اجزای محدود و هم‌هندسه قابل 
پیش‌بینی است؛ زیرا در حل مسائل از شتاب. سرعت 
و جابه جایی گام زمانی قبل استفاده می‌شود. لذا در 
بررسی انجام‌شده به روش اجزای محدود با 
شبکه‌بندی‌های مختلف دیده شد که با افزایش تعداد 
گره‌ها و در نتیجه ِ شدن جزهءهاء نمودار تغییر 
مکان- زمان گره‌های مشخص, دیگر بر هم منطبق 
کیت اما زک ای وم ی کل 

تمام کارهای ارائه‌شده در ایین مقاله در حبطه‌ی 
مسأله‌های همگن قرار دارند و به عنوان کارهای 


پیشنهادی می‌توان این بررسی‌ها را در رابطه با 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


۱ مقایسه روش هم هندسه با اجزای محدود در تحلیل دینامیکی... 


مسأله‌های همسانی عمودی نیز به انجام رساند. مسأله‌های همکن و همسانی عمودی با حضور میرایی 
همچنین با توجه به در اختیار داشتن ماتریس جرم و انجام داد. 
سختی می‌توان کلیه کارهای فوق را در حیطه‌ی 


وس 

وکاصعصهاه ماتص رنلهت ملد رلفمه متامموممع0و1 و.ل۷ ر22116۷5ظ 24 هل متام ۲.1 مقعطمن۳ ...1 
لح ممنمهطع/۱ 0عناموه ما وهمطاه تعنامم رکصه‌صمصه: عمج مه بتاعمعم آمهعه رفظلا۱ 
(2005) ,4135-4195 .00 ,194 ,۷۵1 ,"۳۵1066۲158 

نصماممماصا منتمناافل‌تنا متتامجهممع 190 و۲۰ رعصقط2 200 نک 1.۲ رقمطمنا! پب۷۰۸ ,ملق .۷ بو7ع2211ظ ...2 
(2008) 3-37 .00 ,43 ,۷۵۱ روعلصهط۱۷۱۵0 لمصمتاهانام‌مت و قصمتاهانام‌مصمی 4ص رفصطاتتمعل2 رنوتمعطا 

,"1106 طا تعتفصه احمط ۵۶ فلوررلفمه متتاممصممع 19۵ ر.م) رعصقط2 24 .لا ره ۷۷ منک و۷۱0 ون. ۷۷ م2020 .3 
(2012) ,2174-2178 00۰ ,105 ,۷۵۱ ,۱۷۲۵۲6۵11215 20 ومنصهطمع۱۷( ۵0۵1160 

ماه ۵۶ فلوولفمه متتمصممع۵و1 ونگ.ل.۲ مققطفن۳۳ مضه ۷۰ ر22116۷5ظ یب ملاع یم بلاه‌تاهن .4 
,5257-5296 00۰ و195 ۷۵۱ رعطلت6م۲۵۵1 مه ومنصهط۱60 1160ممرظ 1 ف0مطام۱ امن و قط۱۷۱0۲200 
2006(۰) 

۶ متام را قصمتهط۱1 اهامای که رنه مطا رمک طمهمتومره فلو رفظ متتامجممع150 م۸" ربق رتلقع ...5 
(2006) ,1-30 00۰ و10 ۷۵۰ رعصلتععصنم در معلمتمطانه۴ 

1 _ ابوالبشری م.ح. حسنیء ب. حسن‌زاده» ع. ژیانی عیدگاهی, ا. و ظریف مقدم. ن.» «تحلیل ارتعاشات آزاد مسائل تنش صفحه‌ای با 

استفاده از روش ایزوژئومتریک» ‏ اولین کنفرانس بین‌المللی آکوستیک و ارتعاشات. دانشگاه صنعتی امیرکبیر» تهران ایران» (۱۳۹۰). 

10 همم ۵۶ ممتامی آهمنتاممم۵ ۰150۵ رب۸.٩‏ رتامف12721 مضه ۱۱.2۰ مرصصمل20طع۱۷۱۵ ربظ رتصهععها ...7 
(2009) ,579-592 00۰ ,6 .310 و10 ,۷۵۱ و(عصافبامط مه مصتآنباها) مصلهمصلم مه اه ۵۶ امصیامز صواو۸ 

۸ حسنی, ب. و ظریف مقدم. ن.«توسعه روش عددی جدید برای حل معادلات دیفرانسیل متداول با استفاده از توابع پایه اسپلاین‌ها» 

گزارش فنی شماره ۱۰۱۵ دانشگاه صنعتی شاهرود ایران» (۱۳۸۸). 

1993(۰) ۷0۲1 ۱۱۵ رألن-۵02۲] ۷ رفمتناامناو ۵۶ ومتصصهطونا ور رهم۳۵۵2 2 ۱:۲۰ رطعتامان ,9 

51206 نامع و[ ۱ و مصقللهمطوه۱۷]۵ مج ,9.۷ رتاف2/21] ب مل22ص2طک ر.ظ متصعععع۳۱ .10 
مج ومد مه وومتمممن ۲۷۷۵۲۱۵ . طاق .. وتقصهماطامم. لقصمتعصمصت.. معط ۵ ممتاهمتصتامن 
(2009) مله۳۵۲6۵8۵ رط0حاولن ,1-5 مصیال رطمتاه2نصتاون0 تفصتام معا ات۷۲ 

,60 ۹66000 ,(مونامه‌تمتصهمی) لمناع1 ۷ ط1 وطاورهنع۱۷۲۵۵۵) ۲۱۲۲5( ۲۲۵ ,۷۷۰ ,111167 200 .1 ملع(۳ .11 
٩01۲820-۷712, ۵۷ ۷0۲1, )1997(۰‏ 

مع۷۵ ر.ع 60080 | آممتممن ۱ وتا ما ممتامنل‌مطوا و۲۸ .1.۳ ,۴0805 .12 
2001(۰) و0۲6ط5اهانام صصعحصگیتی 

1993(۰) و1960 ملأ۱۷]۱۵0۲۵۷-۳11 ,۵0( اصعصصعاظ ماتصز؟ عمط ما ممتامیه00ع] و۸ ۲۲ ولا . ,۷ ۳۴۳12 .13 


نشریه مهندسی عمران دانشگاه فردوسی مشهد سال بیست و چهارم. شماره یک, ۱۳٩۱‏ 


